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1 Introduction 



The starting point for most calculations of nonleptonic kaon decays is an effective 
weak lagrangian of the form pj, || 

Cl(\AS\ = 1) = y/2G F V vd VZ,'£C i O i (1) 

i 

which can be derived with the help of the Wilson operator product expansion from 
elementary quark processes, with additional gluon exchanges. In the framework of 
perturbative QCD the coefficients Cj are to be understood as scale and renormaliza- 
tion scheme dependent functions. There exist extensive next-to- leading order (NLO) 
calculations |4], || in the context of kaon decays, among others. These calculations are 
based on the possibility of factorization of short- and long-distance contributions into 
Wilson coefficient functions Cj and mesonic matrix elements of four-quark operators 
Oi, respectively. The latter, however, can presently be obtained only by using non- 
perturbative, i.e. model-dependent, methods, because not only perturbative QCD 
breaks down at scales \i < lGeV, but also the QCD degrees of freedom (quarks and 
gluons) have to be replaced by the mesonic ones. Thus, a fully satisfactory solution 
would include the theoretical understanding of confinement. The only consistent 
approach to this problem may be found in lattice calculations. A discussion of the 
present status has been given in || and will not be repeated here. 

Usually, the results of calculations are displayed with the help of S-factors in the 
form 

T K ^2n = V2G F V ud V* s [Ci{iM)Bi(fMj\ < nn\Oi\K > vac . sat . , (2) 

i 

where the mesonic matrix elements of four-quark operators are approximated by their 
vacuum saturation values, which are real and /i-independent. In principle, factors 
Bi(fi) should be estimated by some higher-order calculations in the long-distance 
regime, for instance, in l/iV c -expansion ]7| in the form 1 + 0(1/ N c ), or from the 
lattice approach. The preliminary stage of these calculations is best characterized by 
the long standing difficulties to explain quantitatively the well-known AI = 1/2 rule. 
Of course, the lack of such calculations for long-distance effects severely restricts the 
predictive power of ([!]), leaving only the possibility for some semi-phenomenological 
treatment ||, [T], with correspondingly large theoretical uncertainties. 

The main aim of the present paper is a further semiphenomenological treatment 
of the long-distance (non-perturbative) aspects of the above lagrangian, especially 
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in view of the actuality of the task to analyze the implications of the measured pa- 
rameter of direct CP violation, e je. The features distinguishing our approach from 
others [fH |9j [Til p], f[2[l are mainly the following: 

- At first, a chiral lagrangian up to 0(p 8 ) is used for deriving mesonic currents and 
densities, from which the matrix elements up to 0(p 6 ) are constructed. 



- At second, according to Weinberg's power counting scheme the calculation 
includes tree level, one- and two-loop diagrams, whereby the renormalization of the 
perturbative (loop) expansion for the matrix elements in question makes use of super- 
propagator regularization (to be discussed below), being connected to the intrinsic 
scale AttFq of the chiral lagrangian, and showing good stability (decreasing higher 
order contributions). 

While the consistency of this approach could be checked phenomenologically in the 
strong interaction (hadronic) sector, in the application to the weak AS = 1 interac- 
tions there remains the above problem of matching the scale- and scheme dependence 
of the Ci, thereby bridging the gap between short and long distance treatments. 

In the present paper, we perform the calculation of matrix elements successively 
for increasing orders in p 2 , displaying also the intermediate results, in order to analyze 
the trend of the successive expansion terms. An obstacle to this procedure is the 
proliferation of structure constants in higher order chiral lagrangians, which have 
to be fixed by experiment. This has been accomplished up to now only to 0(p A ) 



14, [L5L [L6|. As a way out, we invoke another effective model - the Nambu-Jona- 



Lasinio (NJL) model JlTj - whose modifications have been used by several groups 



IS, 191 EOJ to "derive" the chiral lagrangian by bosonization of the fermionic degrees 



of freedom, suitably adapting the free parameters to reproduce those of the chiral 
lagrangian. In this framework, the structure constants of higher order lagrangians 
can be calculated, and they are well comparable in the 0(p A ) case with the empirical 
ones. Therefore it seems justified, to estimate effects of orders beyond 0(p A ) by 
taking NJL-derived structure constants. 

If we compare the predicted amplitudes (g) with experiment, it turns out that 
even after replacing < 7C7t\Oi\K > vac . sa t. by higher order matrix elements there are 
still some correction factors needed, which we call Bi. We restrict ourselves in this 
paper to the display of their ranges and correlations, especially taking into account 
the large value e obtained in the NA31 experiment [2~I] and confirmed recently by the 
KTeV [f22J and NA48 |23[] collaborations. As our approach to the renormalization 
of chiral perturbation theory involves no arbitrary cut-off or scale other than Fq, 
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the bare tt/K decay constant, there is also no other possibility to match the scale 
dependence of the Wilson coefficients, except that due to the renormalization of Fq 
(and the other bare parameters of the effective lagrangian). This is at least partly 
included in our approach, as we redefine the bare coupling F for each order to reach 
agreement with tt/K — > \w decay As mentioned above, this procedure is stable and 
consistent, i.e. does not lead to large higher order corrections or large renormalization 
effects. 

In section 2 we repeat all relevant definitions taken from our earlier work. Section 
3 discusses the higher order structure constants used in the calculation of K — > tttt 
amplitudes, the latter being sketched in section 4. The last two sections give our 
results and conclusions. 



2 Lagrangians and currents 

In the present paper we use the operators 0; in the representation given in [f2|, [24}]: 

2 = ul^^ul diYsi + dhl^u L ulYs l + 2d L 'y fl d L d L ^s L + 2s L ^^s L d L ^s L , 

03 = ul'J^ul d L ^s L + Jl7 m m l ul^sl + 2d L ^d L dj/fs L - 3s L j^s L d L ^s L , 

04 = u L ^^u L d L Ys L + dLl^UL u l Ysl - dLl^d L d L ^s L , 

C>5 = d L ^X a c s L l Y qr^Kqr) > & = d Lltl s L [ Y QRl^QR 

\q=u,d,s J \q=u,d,s 

7 = 6d L ^s L l Y Qr^QQr) , 8 = 6d L ^\ a c s L i Y QR-fKQQR 

\q=u,d,s J \q=u,d,s 

where q^n = \ (1 =F 75)?; K are the generators of the SU(N C ) color group; Q is the 
matrix of electric quark charges. The operators 05^ containing right-handed currents 
are generated by gluonic penguin diagrams and the analogous operators 07,8 arise 
from electromagnetic penguin diagrams. The operators 01,2,3,5,6 an d 04 describe 
the transitions with AI = 1/2 and A J = 3/2, respectively, while the operators 07,g 
contribute to the transition with both A J = 1/2 and A J = 3/2. 

The Wilson coefficients Cj of the effective weak lagrangian (JJ) with four-quark 
operators 0j are connected with the Wilson coefficients Cj corresponding to the basis 
of four-quark operators Qi given in Refs. [1], [|, by the following linear relations: 

Ci = Ci — C2 + C3 — C4 + Cg — C10 , C*2 = r(ci + C2 — Cg — C10) + C3 + C4 , 

5 
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c 3 

C*5 — c 6 



^C*4 = ^(jj( C l + C 2) + C 9 + Cio 



(3) 



The bosonized version of the effective Lagrangian (|l|) can be expressed in the 
form |25|: 

1 



C 



G F 



(4m — *4/J(^z^+ iJln) ~ (4/j+y^4^)(4/j + iJifi) 

+(6 + 5 6)^|^(JS M + iJ£„) + -^£ 3 4^(4 M + »^ M ) 



+^4 
"4^5 



'Lp, bo 'LfiJK° 'Lfj,~ r bu Ln! "T ^ u LfiK u Lfi^ lu 

(jR-iJ 2 R )(Ji+iJt) ~ (4-^4- /f 4)(4+*4) 



'-(J R +iJ 7 R )(V2Jl-J\ 



-166 



(4-*4)(4+*4) + o(4 



4- 



4)(4+*4) 



+ h.c. 



(4) 



Here G F = \f2G F V ud V* s , J^/r^ an d </£/# are bosonized (V =F A) and (S =F P) 
meson currents and densities, corresponding to the quark currents 57^4(1 =F7 5 )A a g 
and densities g|(l =p 7 5 )A a g, respectively (A a are the generators of the U(3)p flavor 
group); 



6 = Ci [1 - 

£5,8 = C5.8 1 



f 



4, + =U, 



AT 2 



^2,3,4 — ^2,3,4 ^1 "+ 
^6,7 



1 



2N, 



(5) 



where the color factors 1/N C originate from Fierz-transformations of four-quark op- 
erators Oi (see more technical details in p5| ). 

Only the even-intrinsic-parity sector of the chiral strong lagrangian is required 
to describe nonleptonic kaon decays up to and including 0(p 6 ). The meson cur- 
rents/densities J1/ RfJL and J1/ R are obtained from the quark determinant by variation 
over additional external sources associated with corresponding quark currents and 
densities [EH]. From the momentum expansion of the quark determinant to 0(p 2n ) 



one can derive the strong lagrangian for mesons C e f f of the same order and the corre- 
sponding currents and densities Jl/r^ anc ^ ^l/r ^° the order 0{p 2n ~ l ) and 0(p 2n ~ 2 ), 
respectively. For example, from the terms of quark determinant of 0(p 2 ) one obtains 
the following: 



J 



(p 2 ) 

eff ~ 

(P> 
L/j, 



-T tr(L ' ) + ? tr(xf/t + f/xt: 



iF 2 



pj i 



J 



(p°)a 



F 2 

-^mR tr (\ a U) 



(6) 



where U = exp (^r~¥>) , with ip being the pseudoscalar meson matrix, and = 
DJJlfl, Dfl = dpU + (Af^U - UAf) and A^ L = V ll ±A ll are right/left-handed 
combinations of vector and axial-vector fields. Furthermore, F « 90 MeV is the 
bare coupling constant of pion decay, X = diag(Xu, Xdi Xs) = —2m <qq> Fq 2 is the 
meson mass matrix, xl = 0.0114 GeV 2 , xl = 0.025 GeV 2 , xl = 0.47 GeV 2 , m is the 
current quark mass matrix, <qq>= (— 220MeV) 3 is the quark condensate, fn ~ 265 
MeV is an average constituent quark mass, and R =<qq> / (juFq) = —4.96. 
At 0(p 4 ) one gets 



C 



(p 4 ) 

eff 



trL 



2\2 



- L 2 n(-[L fl ,L u }' 2 + 3(7; 



2\2 



-U tr (L 2 ) tr ( X U j + U X ] ) - U tr U 2 ( X tf + U X 



x(P 3 )« 
j(p 2 )a 



+L 8 tT[(x j U) 2 + (xU^ 2 )+H 2 tixx j , 

itr {A a [L 4 ^tr( X f/t + ^t) + Il 5 {L m , ( x [/t + ^t)}] J , 

-roiitr {A a [L^ tr(Lj) + L 5 (LjZ7) - 2L 8 U X j U - H 2X ] } 



(7) 



where Li and if 2 are structure constants introduced by Gasser and Leutwyler |14j| . 

For the sake of brevity, here and in the following expressions for the lagrangian at 
0(p 6 ) we restrict ourselves to the terms which are necessary to calculate the decay 
K — > 2n. At 0(p 6 ) one needs the following terms: 1 : 



C 



(p 6 ) 

eff 



tT\Q 12 (xR^(D^D v U + D v DJJ)rfV 

+X ] L»U(D I1 D U U ] + D V D^)UR U ] 
+Qiz[x(D ll D u rfL»L» + RVR^UD^Dvrf) 
+X ] (D^D V UR»R U + UWD^DvU) 



1 The rather lengthy full expression for the bosonized effective lagrangian at 0(p e ) was presented 
in Refe.Ecfi. 
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+Qu W^DpDJjWffrt + D^D l/ U j D»D v UU^ 
+X j (UD fl D l/ U j D^D u U + DpDvUDVtfu) 

+Qu>x%xR' i + Qiz(x ] xR^ + xx ] L^) 
+Qw[(xU j ) 3 + (x j u) 3 } + Q 20 (u j xx j x + u x ] xx' 



where Qi are structure constants introduced in POfl , whereas = U^D^U. The 
corresponding terms of (V =F A) and (S =F P) bosonized meson currents are given by 



J, 



(P 5 )a 



i-tr|A ra [- 2Q 14 [(U X j + X U ] )D^D U U U ] L» + D^D U U(U^ X U ] + X ] )L U 
- UD U ((U^ X + X ] U)D U D^ + D U D^(U X ] + xU ] )) 
+ L"U((U^x + X j U)D^D v U j + D.DMHU^ + xU ] )) 
+ D V ((U X * + XU ] )D V D,U + D U D,U(U^ X + X ] U))u\ 

+2Q l5 {U X ] L, X U ] + XU ] L,U X ] ) + 2Q lG {{U X ] xU\ L,} + { XX \ L,}) 

+2Q 17 ({(Ux^,L lx } + {(xtf) 2 ,L„}) 



and 



jf> mRti[\ a [Q l2 D l U{D^D v }U ] UR U 
+Q 13 (L v L»D fi D I/ U + D„D V U ■ R»R U ) 
+Q li {UD v D tl U ] D u D^U + DvDpUffVrf ■ U) 
+Q 15 L^ X R, + Qw (xRl + Llx) + Q 17 (U X j URl + tift^U) 
+2Q 18 L^L.U + Q l9 (U X ] ?U + Q 20 (xU 1 X + XX ] U + U X ] X)] } ■ 

We do not show explicitly the terms of the effective action at 0(p 8 ) generating the 
scalar current jjf ■* which is necessary for the full calculation of the tree-level matrix 
elements at 0(p 6 ) for the penguin operators, since the corresponding contributions 
turn out to be negligibly small. 



3 Structure constants 

For numerical estimates of K — > 2n amplitudes and e je we will need the values 
of the structure constants Lj and Qi which were introduced in the effective chiral 
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lagrangians at 0(p A ) and 0(p 6 ), respectively. The current experimental status of 
the effective chiral lagrangian at 0(p A ) has been discussed within ChPT in some 
detail in Hl6| . For the 0{p A ) lagrangian ([/]) all structure constants Lj are at present 



determined phenomenologically as measurable values L[ depending on the renormal- 
ization scale jl. The best values of the parameters Lj quoted at a p-meson mass scale 
and the sources of the experimental information used are listed in table 1. The scale 
dependence of the measurable coefficients L\ is determined by relation 

L r M = L\{m) + 7^ ln - , (9) 
(4tt) 2 /i 2 

where the coefficients r« are also given in table 1. 

In the context of the scale dependence of the structure coefficients Li, we have to 
note that in our approach the UV divergences resulting from meson loops at 0(p A ) 
and 0{jP) were separated by using the superpropagator (SP) regularization method 



[ 26(1 which particularly well suits the treatment of loops in nonlinear chiral theories. 
The result is related to the dimensional regularization technique though some the 
difference lies in the scale parameter jl which is no longer arbitrary but fixed by the 
inherent scale of the chiral theory jl = AttF q « 1 GeV, and the UV divergences have 
to be replaced by a finite term using the following substitution: 

(C - l/e) -> C SP = -1 + 4C + (3ir, 

where C = 0.577 is Euler's constant, e = (4 — D)/2, and (3 is an arbitrary constant 
introduced by the Sommerfeld- Watson integral representation of the superpropagator 
based on unitarity. 

The phenomenological analysis of the so-called Skyrme and non-Skyrme struc- 
tures in the effective chiral lagrangian at 0(p A ) was earlier carried out in [E8| by 



using the direct SP-calculations of meson loops for 7T7r-scattering amplitudes. After 
reformulating this analysis in terms of the structure coefficients Lj, the values 

L x = (0.6 ±0.2) • HT 3 , L 2 = (1.6 ±0.3) ■ 10" 3 , L 3 = (-3.5 ± 0.6) ■ 10" 3 (10) 

were obtained from the experimental data on 7T7r-scattering lengths. In the same 
way, taking into account the tadpole loops, the splitting of the decay constants F n 
and Fk was used at 0(p A ) to fix Csp ~ 3.0 and L 5 = (1.6 ± 0.3) ■ 10 -3 . The latter 
value as well as the values (|To"D are in a good agreement with the corresponding ones 
given in table 1. This fact indicates that the choice jl = m p for the renormalization 
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scale of ChPT proves to be consistent with the internal scale of SP-regularization. 
Therefore, we use the values of Li given in table 1 for further phenomenological 
analysis. 

The structure constants Qi of the 0(p 6 ) lagrangian (§) are still not defined from 
experiment. Therefore we need some theoretical model to estimate their values. 
Both the structure constants Li and Qi can be obtained from the modulus of the 
logarithm of the quark determinant of the NJL-type model [|17| which explicitly 
contains, apart from the pseudoscalar Goldstone bosons, also scalar, vector and axial- 
vector resonances as dynamic degrees of freedom. However, in order to avoid double 
counting in calculating pseudoscalar meson amplitudes when taking into account 
resonance degrees of freedom, one has to integrate out (reduce) these resonances in 
the generating functional of the bosonization approach. As a consequence of this 
procedure, the structure coefficients of pseudoscalar low-energy interactions will be 
quite strongly modified. In this way one effectively takes into account resonance- 
exchange contributions J[9|, ^7], ^(J . 

Without reduction of resonance degrees of freedom, the structure constants Lj = 
N c /(16ir 2 ) ■ li, and Qi = N c /(327r 2 m 2 ) ■ qi are fixed from the bosonization of an 
NJL-type model as 

11 1 

^i = 7^^2 = 7TT ' ^ = _ c> ^4 = 0, l 5 = xy-x, 
Z 24 o 

k = ^xy -x y - — , 

and 

1 1/1 \ x 

Ql2 = i Ql3 — X + C) , 014 = — , 

2 . . /I \ 1 4 o x. , / 1 

<7l5 = —Xll — X) — 2x C, Oik = 1 X H (1 — 4x1— 2 £ 

y 5 3 V ; V3 / ' y 6 120 3 6 V ; V 6 

1 x „ A \ ( 1 \ 4 2 f 1 \ 
an— 1 — (1 — 4a;)— xH — c, oia = -x + x c, 

y 120 6 V ; V 6/ y 3 V6 / ' 

1 2 
919 = - TTTjj -x 2 + -x 3 + x(l + 2xy)c , 

qio = 7^5 + X* + 2(1 - 2y)x 3 - x(l + 2zj/)c , 

where x = -mF$/(2<qq>) = 0.1, y = 4tt 2 F 2 / \N c m 2 ) = 1.5 and c = 1 - l/(6y). 
After reduction of the resonances, the structure coefficients get the form 



jred 1 ired ^ 

^ ~2~' 2 "12 



, v + 2(^-l)(ly(^-l)-Zj 



S 



ired 



ired 
M 



Z» A +3(Z 



7 1 



0, z 



red 



(S 



.4 i 



red 



16' 



^1 



7 2 



and 



Qlf = 




q[? = - 


64 r 


-Wif = 


- 1 

~48 A 


-n red — 
3 Qis ~ 


- 1 Z 


192 



24' 



',4 ) 



7 2 



4- 



2Zi 



^(3y - 2) 



6(1 + 4(1 

5-12(1 



1- 4(l + 16(1- Z 

(i-y) 2 



1-2/ 



where y = 47r 2 F 2 /(Z^iV c m 2 ) = 2.4, and = 0.62 is the n — A\ mixing factor. 

In table 1 we also present the predictions of the NJL model for the structure 
coefficients L, which after reduction of meson resonances turn out to be in a good 
agreement with phenomenology. This fact indicates that the NJL-model is a reason- 
able low-energy approximation for the effective four-quark interaction, generating 
a realistic effective meson lagrangian. Therefore we also use it to fix the values of 
the structure constants Qi for numerical estimates of the contributions of the 0{p & ) 
lagrangian (|8|). 



4 Amplitudes of K — > 2n decays 

Using isospin relations, the K — ► 2n decay amplitudes can be parameterized as 

[2 1 [2 2 

Tko-^w+tt- = y g A + -j= A 2 , T K o s ^ 0lT o = y - A - -j= A 2 . 

The isotopic amplitudes A 20 determine the K — > 2n transitions into states with 
isospin 7 = 2,0, respectively: 

A 2 = a 2 e i5 \ A = a e iS \ 

where #2,0 are the phases of 7T7r-scattering. It is well known that direct CP violation 
results in an additional (small) relative phase between a 2 and a$. Let us next in- 
troduce the contributions of the four-quark operators Oi to the isotopic amplitudes 
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A^ by the relations 

8 



Aj = TjAj, Aj^-i^dA?, (11) 

i=l 



where T 2 = y/2T Q = ^G F F (m 2 K - m 2 ). 

At 0(p 2 ), corresponding to the soft-pion limit, for the nonzero tree-level ampli- 
tudes Aj we obtain the following expressions: 



4 1} = -4 2 ' 3) = -4 4) = -1 , 4 7) = -4 7) = 2 , 4 5) = -32(^) L 



(8) _ 16(i?m) 2 f i 2 



^(8)_ 8(/?m) 2 f 2 
- m 2 I F 2 



o 

+ (L 5 - 4L 8 )( X 2 + 3^ + 2x1) + 2L 5 ml 
+ (L 5 - 4L 8 )( X 2 + 3^ + 2x1) + a^m 2 , 



(12) 



The L$ and i?2 contributions in the penguin operators 0§ : % also have a tadpole 
contribution from K — ► (vacuum), included through strong rescattering, — > 7171".^" 
with K — > (vacuum). At 0(p 2 ), in case of the penguin operator 5 , the L 8 and 
i?2 contributions to the direct matrix element from K — > 27r vertices, are fully 
canceled by the tadpole diagrams 0. This is due to the possibility of absorbing the 
tadpole contribution into a redefinition of the K — > 27r vertex if all particles are 
on mass shell. Moreover, such a cancelation is expected at all orders of K — > 2n 
amplitudes including loop diagrams due to general counter term arguments given 
fL5| . According to these arguments the structure constant H 2 is not directly 



in 



measurable and does not occur in the amplitudes of physical processes. 

Some interesting observations on the difference of the momentum behavior of 
penguin and non-penguin operators can be drawn from power-counting arguments. 
According to Eq. ([]) the leading contributions to the vector currents and scalar den- 
sities are of 0(p r ) and 0(p°), respectively. Since in our approach the non-penguin 
operators are constructed out of the products of (V — A)-currents J]- , while the pen- 
guin operators are products of (S — P)-densities J£, the lowest-order contributions of 
non-penguin and penguin operators are of 0(p 2 ) and 0(p°), respectively. However, 
due to the well-known cancelation of the contribution of the gluonic penguin operator 



2 We thank W.A. Bardeen and A.J. Buras for drawing our attention to this point. 
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(9 5 at lowest order [p9|| , the leading gluonic penguin as well as non-penguin contribu- 
tions start from 0(p 2 ) []. Consequently, in order to derive the (V — A)-currents which 
contribute to the non-penguin transition operators at leading order, it is sufficient to 
use the terms of the quark determinant to 0(p 2 ) only. At the same time the terms of 
the quark determinant to 0(p A ) have to be kept for calculating the penguin contribu- 
tion at 0(p 2 ), since it arises from the combination of (S — P) -densities from Eqs. (^) 
and (0), which are of 0(p°) and 0(p 2 ), respectively In this subtle way a difference in 
momentum behavior is revealed between matrix elements for these two types of weak 
transition operators; it manifests itself more drastically in higher-order lagrangians 
and currents. This fact makes penguins especially sensitive to higher order effects. 

Our calculations involve Born and one- and two-loop meson diagrams and take 
into account isotopic symmetry breaking (7r°— 77— 77' mixing). The use of a specialized 
analytical computation package based on REDUCE 0] to calculate amplitudes and 
loop integration makes it possible to evaluate a large number of loop diagrams arising 
for different charge channels. The main problem in the calculation of amplitudes at 
0(p 6 ) is the evaluation of two- loop diagrams. A part of them is shown schematically 
in figure 1 (we do not show rather trivial diagrams with tadpole loops). The diagrams 
of figure la were calculated analytically, because the integration in every loop can 
be performed independently when using the superpropagator regularization. The 
two-loop diagrams of figure lb,c,d cannot be calculated analytically, but they can 
be estimated numerically through a dispersion-relation approach in the same way 
as it was already done in [[31]] for the so-called "box" and "acnode" diagrams. Such 
numerical estimates have shown that the contributions of diagrams of lb,c,d do not 
exceed 2% and can be neglected. 

Table 2 presents the modification of the amplitudes A.j when including succes- 
sively the higher order corrections at 0(p A ) and 0(p 6 ). In our numerical estimates 
the Born contribution at 0(p A ) and the one- loop contribution at 0(p 6 ) were calcu- 
lated for central values of the phenomeno logical parameters Lj from table 1. The 
Born contribution at 0(p 6 ) has been estimated for values of structure constants Qi 
fixed from the bosonization of the NJL-model with reduction of meson resonances. 
Table 2 shows that the Born contribution at 0(p 6 ) is very small as compared to loop 
contributions and does not play an essential role in our further analysis of decay 



3 There is no cancellation of the contribution of the electromagnetic penguin operator Og at the 
est 
0(p°). 



lowest order and the first terms in the expressions (O) for correspond to the contributions at 
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amplitudes and e /e. 

A strong indication that the development to higher orders is physically sensitive 
is given by the behaviour of phases: the strong interaction phases #2,0 arise first at 
0(p 4 ), but for the quantitative description of the phases it is necessary to go beyond 
0(p 4 ). At 0(p 4 ), for the 7T7r-scattering phase shifts and their difference A = 5q — 62, 
we have obtained the values of 5q ~ 22°, 62 ~ —13°, A m 35° which are in agreement 
with |||. At 0(p 6 ), however, we obtained S « 35°, S 2 ~ -9°, A « 44°, in a better 
agreement with the experimental value A exp = (48 ± 4)° [B3 



5 Phenomenological results 

In our approach the parameters £j in Eq. (|TTD are treated as phenomenological (/i- 
independent) parameters to be fixed from the experimental data. They can be related 
to the //-dependent QCD predicted by using some //-dependent Sj-factor defined 
as 

The factors Bi(fi) can be related to the factors defined in (0) by obvious 

relations. Table 3 shows the QCD predictions for the coefficients = CI^CaO + 
t^ v \h) which correspond to the Wilson coefficients 

V V* 

Ci(ji)=Zi(p)+Tyi(p), t = - * l l , 

V udV us 

from the table XVIII of Ref. |1| calculated numerically from perturbative QCD 
at [L = 1 GeV for m t = 170 GeV in leading (LO) and next-to-leading orders in 
"naive dimensional regularization" (NDR) and t-Hooft-Veltman (HV) regularization 
schemes. The numerical values of the QCD scale A^- given in table 3 correspond to 
aj^(M z ) = 0.119 ± 0.003. ^ and Q v ' were obtained from Zi and respectively, 
using the Eqs. (|) and @. 

As we cannot calculate the factors -B«(/i) theoretically, they can be fixed only 
from data in the spirit of the semi-phenomenological approach |], |], §]. Table 2 
shows that the amplitudes of K — » 2ir decays are dominated by the contribution of 
the operators (9; with i = 1,2, 3, 4, 5, 8. Moreover, in case of the operators ©1,2,3, the 
first term in the combination (— £1+ £2+ £3) dominates in the effective weak meson 
lagrangian (f|). Thus, the isotopic amplitudes can be given after restriction to the 
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dominating contributions of four-quark operators as 



+& y \»)BM 4 5) + d^V)^) 4 8) ] , (is) 

and the relation ([11]) to the measurable amplitudes may be modified to 

A!={aT + raf)e lSl . 

At least two factors B\ and P 4 can be estimated from the experimental values Aq XP ~ 
10.9 and A^° p ~ 0.347 (for fixed P5 and Pg) while the other two (penguin) factors 
_E?5 and B$ should be fixed from other data. The factors Pi, B4, P5 and B$ are the 
analogs of the bag factors B^ 2 \ b[ 3 ^ 2 \ B^ 2 ^ and B^ 2 \ respectively, introduced 
in0. 

The parameter e of direct CP- violation in K — > 2tt decays can be expressed by 
the formulae 

u Ima 0/1 _ iU/o +s ^s^ Rea 2 „ llma 2 



i« i A z ^)(,.\u /i( 4 ) 



e = — - 



a/2 Re a Re a ' u Im a ' 

and the ratio e je be estimated as (recall that, experimentally, s /e ~ Ree'/^j arge ~ 
arge') 

(3/) 

1 = Im A, (P -P 2 ), Pj= -^ , " , ^y, (14) 

with ImAt = ImV^*Vtd = \V u b\ \ V c b\sin5 = r]\V us \\V c b\ 2 in the standard and the Wolfen- 
stein parameterizations of the CKM matrix. 

Table 4 gives the estimates of e /e from a semi-phenomenological approach ob- 
tained after fixing the correction factors Pi and P4 for isotopic amplitudes in the 
representation (|T3|) by experimental (CP-conserving) data on Re .4.0,2, and setting 
P5 = B 8 = 1. We have used the matrix elements of the operators Oi displayed in 
table 2 (for central values of phenomenological structure coefficients Lj given in table 
1), and the theoretical values from table 3. The values (e /e) m in an d {s /e) ma x 
correspond to the interval for Im A t obtained from the phenomenological analysis of 
indirect CP violation in K — > 2tt decay and P° — P° mixing |l|, ||: 

0.86 • 10" 4 < Im \ t < 1.71 ■ 10~ 4 . (15) 

Table 4 demonstrates the modification of the semi-phenomenological estimates of 
the parameters Pi, P4 and [e je) rnax after successive inclusion of the corrections at 
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0(p 4 ) and 0(p 6 ). Most important are the corrections at 0(p 4 ). The peculiarity of 
the results at 0(p 2 ) lies in the observation that all estimates of e je lead to negative 
values. This is related to the fact that in the case corresponding to table 4a the 
contribution of gluonic penguins to AI = 1/2 transitions appears to be suppressed, 
leading, after the interplay between gluonic and electromagnetic penguins, to the 
relation P$ < P2 for the two competing terms in (|T4]) . Generally speaking, AI = 1/2 
transitions loose importance compared to AI = 3/2 when estimating e je. The 
situation already changes after inclusion of the correction at 0(p 4 ), due to relative 
enhancement of the matrix elements for the operator (9 5 (see table 4b). Taking into 
account the dependence of the Wilson coefficients on the renormalization scheme, 
after including the corrections at 0(p 4 ) and 0(p 6 ) we obtained the following upper 
and lower bounds for e je (see table 4c): 

- 3.2 • KT 4 < £/£ < 3.3 • 1(T 4 , (16) 

where the range characterizes the uncertainty from short-distance physics. 

Our calculations have shown that especially the penguin matrix elements are 
most sensitive to various refinements: higher-order derivative terms in chiral la- 
grangians, the reduction of meson resonances, 7T° — 77 — 77 mixing, and meson loop 
corrections. It should be added that the modification of penguin matrix elements, 
discussed in this note, is much more important for gluonic than for electromagnetic 
penguin transitions. This is obvious from the observation that the latter at the low- 
est order contain terms of 0(p°) which remain unchanged when taking into account 
the additional terms derived from the effective lagrangian at 0(p A ). 

We give some results concerning the dependence of the above semi-phenomenological 
estimates for e je on the choice of the penguin correction factors and B% (figure 
2) and on the values of the structure constants Li (figure 3). In figure 3 we show 
the dependencies of e je on the coefficients L 4 , L 5 and L 8 only, to demonstrate the 
appreciable sensitivity to the variation of these parameters within their phenomeno- 
logical bounds given in table 1. It is caused by the fact, that the coefficients L 4 , 
L5 and L§ appear in penguin contributions to the K — > 2-k amplitudes already at 
0{p 2 ) (see (|T2])) while all other structure coefficients given in table 1 appear in the 
amplitudes of higher orders. 

To study the upper and lower bounds for e je corresponding to the variation of 
the parameters Li and Im At within their phenomenological bounds, we have used the 
so called "scanning" and "Gaussian" methods || . In the first case the parameters Li 



14 



were scanned independently within the intervals defined by their central values and 
errors given in table 1 while the parameter Im X t was scanned within phenomenologi- 
cal bounds ([151) . In the second case we calculated the probability density distribution 
for e je obtained by using Gaussian distributions for the parameters Lj with errors 
given in table 1 [] while for the parameter Im A t using the result obtained in || : 

ImA t = (1.29 ±0.22) • 10~ 4 . 

In tables 5 and 6 the "scanning" and "Gaussian" results for e je are given for different 
values of B 5 (B$ = 1). Figure 4 shows typical probability density distributions for 
£ je obtained in the Gaussian case. 

Our results demonstrate that even after taking into account all uncertainties 
related to both phenomenological input parameters and renormalization scheme de- 
pendence, it is still rather problematic to explain theoretically with B$ = B$ = 1 the 
value of the direct CP- violation parameter Re{e je) = (23.0±6.5) x 10~ 4 measured in 



the experiment NA31 at CERN pi ]. The rather high level of direct CP- violation ob- 
served in this experiment was confirmed by recent measurements of KTeV at FNAL 
@, (28.0 ± 4.1) x 10~ 4 , and NA48 at CERN JEJ, (18.5 ± 7.3) x 10~ 4 . Taking into 
account the result of the experiment E731 at FNAL f34]], (7.4 ±5.9) x 10~ 4 , the world 
averaged value is estimated as 

Re(e'/e) = (21.2 ± 2.8) x 10" 4 . (17) 

Finally, we give some results concerning the factor B 5 required to describe the 
experimental value ([lj]) within our semi-phenomenological approach. In figures 5, 
6 and 7 we show the probability density distribution for factors B\, B^ and B 5 , 
respectively, obtained by using Gaussian distributions for the parameters Li, lm\ t 
and e je. The parameters B\, B 4 and B 5 were defined from the experimental values 
of the isotopic K — > 2n amplitudes Aq, and the ratio e je with B 8 = 1. The 
dispersion of these parameter values in figures 5, 6 and 7 is caused mainly by the 
uncertainties of Li and ImAt, while the experimental error of e je is much less in- 
fluence. Figure 7 demonstrates the necessity for a rather large factor B§. It should 
be emphasized, that for even larger values of B 5 , the contribution of nonpenguin 
operators to the AI = 1/2 amplitude are still dominating (see figure 8). In figure 9 
the probability density plots show the correlations between parameters Bi, B 4 and 

4 With exception of L4 which is not determined experimentally and therefore taking uniform 
distribution inside "theoretical" limits —0.8 • 10~ 3 < La < 0.2 • 10~ 3 . 
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B 5 . The correlations between pairs of parameters Bi, B4 and B 5 , B 4 are caused by 
the isotopic symmetry breaking related with tt° — 77 — r\ mixing. The first plot in 
figure 9 demonstrates the strong correlation between B\ and B4. Due to relatively 
small contributions of penguin operators to isotopic amplitudes of K — > tttt decays 
there are no visible correlations between B±, B 5 and B 4 , B 5 . Figure 10 shows the 
correlations between B$ and B$ calculated for central values of the phenomenologi- 
cal constants Li and ImAt with Re(e /e) = 21.2 x 10~ 4 used as experimental input. 
From this figure one can see that even for B 8 = values of B 5 > 2 are necessary to 
explain the large value of e /e (pp . 



6 Conclusion 

From studying the impact of the recently confirmed large e value on the parameter- 
ization of the hadronic weak lagrangian, including step by step various refinements, 
we have shown the necessity for a rather large gluonic penguin contribution (the 
factor B5 is found well above 1, see figure 7). The large B\ and B§ values may be 
a hint that the long-distance contributions are still not completely understood. An 



analogous conclusion has been drawn in [35], where also possible effects from physics 



beyond the Standard Model are discussed. From the phenomenological point of view, 
there is no difficulty in taking (^) as a bona-fide weak current-current lagrangian with 
coupling constants £j to be fixed experimentally. The problems arise when matching 
these parameters to Wilson coefficients derived from perturbative QCD, which is, of 
course, a necessary requirement. It should be remarked that in our approach there is 
also no convincing argument for the large correction factor B\ (due to the AI = 1/2 
rule); but then we may ask, why B\ and B§ should behave differently: as can be 
seen from table 2, the relative changes of the respective matrix elements in going to 
higher powers of p 2 do not differ very much. 

In this context, one should note recent progress in the estimates of the 5-factors 
with a matching procedure based on higher-order calculations in the long-distance 
regime within the l/iV c -expansion. An essential enhancement of the bag factor for 
the gluonic penguin operator by the 1/N C corrections at next-to-leading order in 
the chiral expansion has been observed in fl36"| , where the value B^ 2 ^ = 1.6 ± 0.1 



has been obtained. The similar value, B^ 2 ^ = 1.6 ± 0.3, arising from 0(p 4 ) chiral 
loop corrections, was obtained in [|Tl], |TT|] within the semiphenomenological chiral 
quark model with values of the quark and gluon condensates fixed by reproducing 
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the AI = 1/2 rule. 

Since our results are very sensitive to the relative contribution of the gluonic pen- 
guin operator, the question of its phenomenological separation in K — > 2n decays 
becomes critical, in the context of the AI = 1/2 rule as well as for a very important 
problem of direct CP-violation. CP-conserving K — > 2ir data alone are clearly not 
sufficient for such a separation. It could be accomplished, on the other hand, when 
taking into account Dalitz-plot data for K — > 3n as well as differential distributions 
for radiative decays K — > 2Trj, K — > n2j described by the same lagrangian (p. As 
emphasized above, the reason for this possibility is found in the difference in momen- 
tum power counting behavior between penguin and non-penguin matrix elements, 
which appears in higher orders of the chiral theory, when calculating various param- 
eters of differential distributions, for instance, slope parameters of the Dalitz-plot 
for K — > 3tt. A substantial improvement in the accuracy of such experimental data 
(mostly being of older dates) would be very helpful for such a phenomenological 
improvement of the theoretical situation for e /e (see |23 for discussion of this point 
and |37|, |SS[ for some recent measurements). 



The authors gratefully acknowledge fruitful and helpful discussions with W.A. Bar- 
deen, A.J. Buras, J. Gasser, E.A. Paschos and P.H. Soldan. 
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TZ 



d) 

Figure 1. Topology of the main two-loop diagrams at 0(p 6 ) (diagrams with tad- 
pole loops on the external lines and in the vertices are not shown). The external 
lines denote the momenta. The internal lines correspond to various combinations of 
virtual pions and kaons in different charge channels. The filled circle denotes the 
week interaction vertex, the open circle corresponds to strong interaction. 
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Figure 2. £? 5 and -Bg-dependence of s'/ £ calculated for central values of the 
phenomenological constants and ImA t = 1.29 • 10~ 4 with = 435 MeV. The 
Independence is calculated with B 8 — 1 and the Independence - with _B 5 = 1. 
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Figure 3. Dependencies of e je on L4, L5, and L$ calculated for central values 
of Im A t and other L r coefficients with = 435 MeV and B 5 = B$ = 1. 
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Figure 4. Probability density distributions for e je with B 5 = B$ = 1. 
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Figure 5. Probability density distributions for B\ with B 8 = 1. 
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Figure 7. Probability density distributions for B$ with B 8 — 1. 
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Figure 8. Probability density distributions for the relative contribution of pen- 
guin operators to the AI = 1/2 amplitude. 
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Figure 9. Typical probability density plots for correlations between parameters 
B u B 4 and B 5 calculated with B 8 = 1 and = 325 MeV in NDR scheme. For 
various values of A^ and different renormalization schemes such plots look very 
similar. 
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Table 1. Phenomenological and theoretical values of the structure coefficients 
Li (in units 10~ 3 ). 





Phenomenology [ 




NJL model 




Ll(m p ) 


Input 






Without reduction 
of resonances 


After reduction 
of resonances 


1 


0.4 ±0.3 


K e /± and irir - 


-> TTTT 


3/32 


0.79 


0.85 


2 


1.4 ±0.3 


K e i and irir - 


-> TTTT 


3/16 


1.58 


1.70 


3 


-3.5 ± 1.1 


K e4: and tttt - 


-> 7T7r 





-3.17 


-4.30 


4 


-0.3 ±0.5 


1/N C arguments 


1/8 








5 


1.4 ±0.5 


Fk/F^ 




3/8 


0.98 


1.64 


8 


0.9 ±0.3 


m K o — m K + 




5/48 


0.36 


1.12 






baryon mass ratios 









30 



Table 2. Isotopic amplitudes of K — > 2n decays 



Operators 


0i 2 03 04 05 06 a0 7 a0 8 


0(p 2 ) 


Soft pion approximation 


ReA^ 
ReA% ] 


-1.000 1.000 1.000 0.000 -9.623 0.000 0.016 1.458 
0.000 0.000 0.000 1.000 0.000 0.000 -0.016 0.654 


Born diagrams with ttq — r/ — r/' mixing 


ReA% ] 


0.004 -0.021 -0.039 0.020 0.119 0.004 -0.001 -0.002 
-0.004 0.021 0.039 -0.002 -0.119 -0.004 0.000 -0.016 


Sum 

P 2 


ReA% ) 
ReA% ] 


-0.996 0.979 0.961 0.020 -9.504 0.004 0.015 1.456 
-0.004 0.021 0.039 0.998 -0.099 -0.004 -0.016 0.638 


0(P 4 ) 


Born diagrams 


ReA^ 


-0.247 0.249 0.236 0.008 -1.626 0.000 0.004 0.037 
-0.003 0.001 0.015 0.249 -0.059 0.000 -0.004 0.008 


ldoop diagrams 


ReA ( Q } 
lmA% ] 
ReA% ] 
ImA^ 


-0.171 0.171 0.111 0.001 -2.072 0.000 0.001 0.188 
-0.482 0.482 0.482 0.000 -4.572 0.000 0.008 0.344 
0.000 0.000 -0.004 -0.149 0.001 0.000 0.001 -0.006 
0.000 0.000 0.000 -0.213 -0.004 0.000 0.003 -0.049 


Sum 


ReA^ 
ImA^ 
ReA% ] 
ImA^ 


-1.415 1.399 1.307 0.029 -13.202 0.004 0.020 1.871 
-0.482 0.482 0.482 0.000 -4.572 0.000 0.008 0.806 
-0.007 0.022 0.050 1.099 -0.157 -0.004 -0.018 0.593 
0.000 0.000 0.000 -0.213 -0.004 0.000 0.003 -0.151 


0(p 6 ) 


Born diagrams 


ReA^ 
ReA^ 


-0.003 0.005 0.005 0.000 0.012 0.000 0.000 0.001 
-0.001 -0.001 -0.001 0.005 -0.004 0.000 0.000 0.000 


ldoop diagrams 


ReA ( o } 
ImA^ 
ReA^ 
Im„4* } 


-0.106 0.107 0.018 0.002 -0.151 0.000 -0.002 0.016 
-0.229 0.232 0.232 0.000 -1.582 -0.001 0.004 0.063 
0.000 0.001 0.002 -0.097 -0.004 0.000 0.000 0.007 
0.000 0.001 0.001 -0.077 -0.001 0.000 0.001 -0.007 


2doop diagrams 


ReA ( o } 
ImA^ 
ReA^ 
ImA^ 


0.202 -0.202 -0.220 0.000 1.753 0.000 -0.003 -0.075 
-0.169 0.169 0.142 0.000 -1.704 0.000 0.003 0.115 
0.001 -0.001 -0.001 -0.036 0.000 0.000 0.001 -0.007 
0.000 0.000 0.000 0.034 0.000 0.000 -0.001 0.006 


Sum 

p 2 + p 4 + p 6 


ReA% ) 
lmA% } 
ReA? 
ImA^ 


-1.322 1.309 1.111 0.031 -11.588 0.003 0.015 1.664 
-0.880 0.883 0.856 0.000 -7.858 -0.001 0.014 1.181 
-0.007 0.021 0.049 0.971 -0.166 -0.004 -0.018 0.566 
0.000 0.001 0.001 -0.256 -0.003 0.000 0.004 -0.140 
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Table 3. QCD predictions for the parameters &(//) = Q (fJ>) + tQ (fj), calcu- 
lated with Wilson coefficients = Zi(fj) + T?/j(/i) at // = 1 GeV for m t = 170 GeV 







= 325 MeV 


A Ks = 435 MeV 




LO 


NDR 


HV 


LO 


NDR 


HV 


LO 


NDR 


HV 


A*) 


-1.286 


-1.061 


-1.165 


-1.443 


-1.159 


-1.325 


-1.624 


-1.270 


-1.562 


S2 


0.187 


0.195 


0.198 


0.172 


0.176 


0.182 


0.157 


0.150 


0.165 


^3 


0.129 


0.143 


0.137 


0.122 


0.137 


0.130 


0.115 


0.131 


0.121 


& } 
S4 


0.645 


0.714 


0.687 


0.609 


0.684 


0.650 


0.573 


0.654 


0.599 


£^ 

S5 


-0.008 


-0.020 


-0.008 


-0.012 


-0.032 


-0.013 


-0.016 


-0.056 


-0.023 


S6 


0.000 


-0.003 


0.000 


-0.001 


-0.007 


-0.001 


-0.002 


-0.021 


-0.007 


£^ /a 
S>7 / " 


0.002 


0.003 


-0.001 


0.004 


0.008 


0.001 


0.006 


0.015 


0.032 


d z) la 


0.000 


0.002 


0.001 


0.001 


0.004 


0.002 


0.001 


0.009 


0.067 


Ay) 
£i 


0.044 


0.038 


0.048 


0.054 


0.048 


0.053 


0.065 


0.060 


0.069 


S2 


-0.028 


-0.029 


-0.030 


-0.029 


-0.033 


-0.030 


-0.030 


-0.033 


-0.030 




-0.002 


-0.002 


0.001 


-0.002 


-0.002 


-0.002 


-0.002 


-0.002 


-0.002 


el" 


-0.009 


-0.010 


0.004 


-0.008 


-0.009 


-0.009 


-0.008 


-0.009 


-0.008 




-0.081 


-0.076 


-0.067 


-0.109 


-0.111 


-0.092 


-0.143 


-0.173 


-0.132 


,(») 


-0.033 


-0.042 


-0.021 


-0.049 


-0.076 


-0.033 


-0.071 


-0.139 


-0.051 




0.033 


0.004 


0.006 


0.044 


0.013 


0.016 


0.057 


0.027 


0.032 




0.031 


0.028 


0.031 


0.043 


0.041 


0.045 


0.058 


0.061 


0.067 
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Table 4. Predictions for the parameters of K — > 2tt decays in the semi-phenomenological 
approach (B$ = Bg = 1). The ratio e /e is given in units 10~ 4 . 



a) At 0{p 2 ): 





= 215 MeV 

MS 


= 325 MeV 

MS 


A^ = 435 MeV 

MS 




LO 


NDR 


HV 


LO 


NDR 


HV 


LO 


NDR 


HV 


Bi 
B A 

Po 
P2 

(e 1 s)rnin 

(e/e) max 


6.82 
0.54 
1.76 
2.88 
-0.2 
-0.1 


7.74 
0.48 
1.21 
2.49 
-0.2 
-0.1 


7.29 
0.51 
0.59 
3.69 
-0.5 
-0.3 


6.26 
0.57 
3.17 
4.41 
-0.2 
-0.1 


7.27 
0.50 
2.73 
4.20 
-0.3 
-0.1 


6.65 
0.53 
2.06 
4.33 
-0.4 
-0.2 


5.71 
0.60 
5.02 
6.37 
-0.2 
-0.1 


6.76 
0.52 
5.91 
7.03 
-0.2 
-0.1 


5.83 
0.58 
4.29 
6.45 
-0.4 
-0.2 



b) Up to and including 0(p 4 ): 





A^L = 215 MeV 


AgL = 325 MeV 


AgL = 435 MeV 




LO 


NDR 


HV 


LO 


NDR 


HV 


LO 


NDR 


HV 


Bi 


4.54 


5.13 


4.85 


4.16 


4.79 


4.42 


3.79 


4.41 


3.86 


B 4 


0.48 


0.43 


0.45 


0.51 


0.48 


0.47 


0.54 


0.46 


0.51 


Po 


3.94 


3.25 


2.41 


6.09 


5.68 


4.58 


8.87 


10.46 


7.87 


P2 


3.49 


3.08 


4.20 


5.21 


5.04 


4.99 


7.43 


8.34 


7.34 


/ £)min 


0.4 


0.1 


-3.1 


0.8 


0.6 


-0.7 


1.2 


1.8 


0.5 


(e / s) max 


0.8 


0.3 


-1.5 


1.5 


1.1 


-0.4 


2.5 


3.6 


0.9 



c) Up to and including 0(p 6 ): 





A^ = 215 MeV 

MS 


A M? = 325 MeV 


A^L = 435 MeV 

MS 




LO 


NDR 


HV 


LO 


NDR 


HV 


LO 


NDR 


HV 


Bi 
B A 

Po 
P2 

{e 1 £)min 
(e / £)max 


4.29 
0.53 
3.93 
3.57 
0.3 
0.6 


4.85 
0.48 
3.23 
3.17 
0.1 
0.1 


4.58 
0.50 
2.40 
4.29 
-3.2 
-1.6 


3.93 
0.56 
6.08 
5.31 
0.7 
1.3 


4.53 
0.50 
5.66 
5.15 
0.4 
0.9 


4.17 
0.53 
4.56 
5.07 
-0.4 
-0.8 


3.57 
0.60 
8.87 
7.55 
1.1 
2.2 


4.17 
0.51 
10.44 
8.51 
1.7 
3.3 


3.64 

0.57 

7.87 

7.43 

0.4 

0.7 
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Table 5. Upper and low bounds for e je (in units 10 4 ) for different values of B§ 
= 1) obtained by the scanning method. 



B 5 


MS' 


LO 


NDR 


HV 




(MeV) 


min 


max 


min 


max 


min 


max 




215 


-3.8 


8.5 


-4.0 


7.5 


-7.2 


3.3 


1.0 


325 


-4.5 


11.9 


-5.1 


11.5 


-6.6 


8.1 




435 


-5.5 


16.1 


-5.9 


19.6 


-6.6 


13.8 




215 


-2.4 


16.3 


-2.7 


14.7 


-5.6 


9.7 


1.5 


325 


-2.6 


22.2 


-3.1 


21.9 


-4.1 


16.9 




435 


-2.9 


29.6 


-2.9 


35.6 


-4.2 


26.7 




215 


-0.9 


24.0 


-1.3 


22.0 


-4.4 


16.2 


2.0 


325 


-0.7 


32.5 


-1.2 


32.3 


-2.5 


25.7 




435 


-0.4 


43.2 


0.1 


51.4 


-1.8 


39.5 



34 



Table 6. Upper and low bounds for e'/e (in units 10~ 4 ) for different values of B$ 
(Bg = 1) obtained by Gaussian method. The limits without brackets correspond to the 
confidence level of 68% while the limits in brackets correspond to the confidence level 95%. 



B 5 


MS 


LO 


NDR 


HV 




(MeV) 


min 


max 


min 


max 


min 


max 




215 


-2.1 


3.0 


-2.4 


2.4 


-4.9 


-0.2 






( -5.1 


7.2) 


(-5.1 


6.4) 


( -8.0 


2.6) 


1.0 


325 


-2.3 


4.3 1 


-2.7 


4.0 


-3.6 


2.1 






( -6.3 


9.9 ) 


( -6.6 


9.7 ) 


( -7.6 


6.5 ) 




435 


-2.6 


5.9 


-2.5 


7.6 


-3.5 


4.7 






( -7.8 


13.4 ) 


(-8.1 


16.6 ) 


( -8.5 


11.7 ) 




215 


-0.5 


7.3 


-0.8 


6.5 


-3.2 


3.1 






( -3.8 


14.4) 


(-3.9 


13.0 ) 


( -6.7 


8.3) 


1.5 


325 


-0.2 


10.1 


-0.5 


9.8 


-1.5 


6.9 






( -4.6 


19.5 ) 


( -4.9 


19.3 ) 


( -5.9 


14.4 ) 




435 


-3.8 


14.4 


-3.9 


13.0 


-6.7 


8.3 






( -5.5 


25.8 ) 


(-5.7 


30.8 ) 


( -6.3 


23.6 ) 




215 


0.8 
( -2.9 


11.8 
21.8 ) 


0.5 
(-3.1 


10.6 
20.0 ) 


-1.8 
( -5.7 


6.6 
14.4 ) 


2.0 


325 


1.5 


16.0 


1.2 


15.8 


0.1 


11.9 






( -3.4 


29.1 ) 


( -3.8 


28.9 ) 


( -4.7 


22.8 ) 




435 


2.4 


21.0 


3.2 


24.9 


1.4 


19.1 






( -4.0 


36.5 ) 


( -4.0 


40.7 ) 


( -4.8 


34.5 ) 



35 



